Doped Stripes in Models for the Cuprates 
Emerging from the One-hole Properties of the Insulator 
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The extended and standard t-J models are computationally studied on ladders and planes, with 
emphasis on the small J/t region. At couplings compatible with photoemission results for undoped 
cuprates, half-doped stripes separating 7r-shifted antiferromagnetic (AF) domains are found, as in 
Tranquada's interpretation of neutron experiments. Our main result is that the elementary stripe 
"building-block" resembles the properties of one hole at small J/t, with robust AF correlations 
across-the-hole induced by the local tendency of the charge to separate from the spin (G. Martins 
et al, Phys. Rev. B60, R3716 (1999)). This suggests that the seed of half-doped stripes already 
exists in the unusual properties of the insulating parent compound. 

PACS numbers; 74.20.-z, 74.20.Mn, 75.25.Dw 



The understanding of high temperature superconduc- 
tors is among the most important open problems in 
strongly correlated electrons. A remarkable development 
in recent years is the accumulation of experimental ev- 
idence compatible with stripe formation in the normal 
state of underdoped cuprates Q]. This includes spin in- 
commensurability (IC) in neutron experiments, results 
believed to be caused by stripes separating 7r-shifted AF 
domains [Q. More recently, it has been shown that the 
stripes are metallic |Q , result compatible with proposals 
of the normal state of x=l/8 cuprates as made out of 
half-doped stripes ^ . Whether stripe formation is bene- 
ficial or detrimental to superconductivity is unclear, but 
it appears that stripes are an important ingredient of the 
normal state that cannot be ignored. 

The theoretical explanation of stripe formation is much 
debated. Early work reported stripes in the t-J (at large 
J/t with 1/r repulsions) and Hubbard (Hartree-Fock) 
models However, these stripes were insulating 

with hole density n/j^l.O, different from the experimental 
n/i~0.5 stripes Recently, considerable progress was 
made when doped stripes were reported by White and 
Scalapino within the standard t-J model [|| (see also Ref. 
[0). In Ref. Q the analysis was performed at couplings 
where two holes form d-wave pairs, and the stripes are 
sometimes described as a condensation of these pairs into 
a stripe domain- wall [^. However, experiments are usu- 
ally interpreted as holes moving freely along site-centered 
stripes In addition, the "extended" t-J model with 
hopping beyond neighboring sites, or the standard t-J 
model with very small J/t, are needed |^,|o| to repro- 
duce the insulator one-hole photoemission (PES) disper- 
sion pd] ]. Thus, understanding metallic stripe formation 
requires further work and searching for stripes in the ex- 
tended t-J model, particularly in regimes without hole 
binding and where the absence of phase separation (PS) 
is not controversial, is important to clarify the driving 
mechanism for these unusual complex structures. 



Building upon previous investigations ||6|,[7| , in this Let- 
ter indications of n?i~0.5 stripes compatible with exper- 
iments |l] are reported in the extended and standard t-J 
models on ladders and square clusters. These stripes do 
not seem composed of hole pairs (although pairs forming 
domain- walls may be present at larger J/t than studied 
here [|j). They also exist in the t-J^ model and us- 
ing classical spins 0, implying that the details of the 
AF spin background are unimportant for its stabiliza- 
tion. Moreover, our most important result is that the 
basic stripe "building-block" exists already in the insu- 
lator where one-hole wave functions have a complex spin 
structure with strong AF correlations across-the-hole, re- 
sembling the stripe spin correlations found here numer- 
ically. These results provide a rationalization for stripe 
formation built upon the one hole properties, in regimes 
where spin and charge are almost separated 

The extended t-J model used here is defined as 

H = J^(Si • Sj - ininj) - ^ tim(cjcm + h.c). 
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where tim is t(=l) for nearest-neighbors (NN), t' for next 
NN, and t" for next next NN sites, and zero otherwise. 
The rest of the notation is standard. The t-Jz model 
is obtained by J^J^ and Si • Sj^ S^iS^j, and t'<0 and 
t">0 are relevant to explain PES data p|-[Tl]]. Here the 
Density Matrix Renormalization Group (DMRG) []6|,p^, 
Lanczos [p^ , and an algorithm using a small fraction of 
the ladder rung-basis (optimized reduced-basis approxi- 
mation, or ORBA |l^) are used. Results are presented in 
(i) the small J/t region with t'=t"=0.0, and (ii) small and 
intermediate J/t with nonzero t' and t" These two 
regions have similar physic s 13 ], and the extra hoppings 
are expected to avoid PS ||5|,p^. Intuitively, t',t" increase 
hole mobility, as reducing J/t does, but also avoid ferro- 
magnetism at small J/t [^. Note also that no coupling 
fine-tuning is needed: the results below appear in a ro- 
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bust region of parameter space. 
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FIG. 1. (a,b) Rung hole density (n(r)) vs rung index r us- 
ing DMRG, with PBC along rungs and OBC along legs, to 
illustrate the n^^CS stripe formation, x is the overall hole 
density, (a) corresponds to a 4x8 cluster with 4 holes. Solid 
(dashed) lines are for the standard t-J (t-Jz) model with J=0.2 
(Jz=0.3), t'=t"=0.0 (inset: same as solid lines but for 4x12 
with 2 holes), (b) Same as (a) but for a 4x12 cluster with 6 
holes, J=0.5, t'=-0.3, and t"=0.0. (c) Hole density at rung 
r, defined now as C(r)=^. (noUi) where the sum is over 
sites belonging to rung r, is an arbitrary site of rung r=0, 
and (no Hi) is the hole density-density ground-state correla- 
tion. The cluster is 4x6 with PBC in both directions, 4 holes, 
J=0.2, t'=-0.35, and t"=0.25 (ORBA with -SxlO*^ states). 
The inset shows ground-state energy vs number of states, (d) 
Distribution of one-hole around a second hole projected at the 
open circle position, for the case in (a) at the indicated rungs 
(running horizontally). Full circles areas are proportional to 
the hole density, (e) S(qa,,7r) vs for the clusters, couplings, 
and densities of (a) and (b). 

In Fig.l, DMRG and ORBA results for 4xN clus- 
ters are shown. In Fig. la the rung density for a 4x8 
(4x12) cluster with 4 (2) holes at small J/t is pre- 
sented. Cylindrical boundary conditions (CBC) are used 
i.e. open boundary conditions (OBC) along legs and pe- 
riodic boundary conditions (PBC) along rungs The 
four holes separate into two groups of two holes, surpris- 
ing result since for a square lattice Jc=0.2 is the critical 
value for hole pair binding in the t-Jz model, and in the 
t-J model Jc is expected to be larger ||l^ . Similar results 
are found in the t-J^ model (Fig. la) and at intermedi- 
ate J/t but with t'7^0, which increases the hole mobility: 
Fig. lb with six holes show the formation of three groups 
of two holes as in Fig. la. This is not spuriously caused 
by the OBC along legs, as shown in Fig.lc with results 
using PBC in both directions. As ORBA starting config- 
uration holes clustered (phase separated) or spread apart 
(free gas) were used, with PBC or CBC, and in both cases 
the results converged to the same "stripe" answer. 

To study the two-hole state internal structure, in 
Fig. Id the density distribution of one hole around the 



other is shown, for one of the 2-hole regions of Fig. la. 
The largest density is at two lattice spacings along the 
rung, and the hole distribution does not resemble a 
tight d-wave bound state ||l^ . Similar conclusions were 
reached for the two holes of Fig.lc. The result is actually 
compatible with the formation of a short site-centered 
stripe where the two holes form a closed loop with den- 
sity 0.5 along a rung |l^. These stripes appear to occupy 
more than one rung in Figs, la-c, and thus they could 
also be labeled as bond-centered However, this effect 
seems to arise from stripe tunneling between neighboring 
rungs, as the one-hole projection suggests (Fig. Id). Sim- 
ilar results regarding half-doped stripe formation were 
also found on 6x6 clusters, as exemplified in Fig. 2a where 
sets of 3 holes form individual iih'^O.b stripes (invariance 
under reflexions was assumed along the legs). Overall 
the results are consistent with Tranquada's description 
of stripes ||l|. They are also consistent with numerical 
reports for the standard t-J model although our in- 
terpretation of the results (below) is different. 




FIG. 2. (a) Rung hole density (n(r)) vs rung index r using 
DMRG on a 6x6 cluster with 6 holes, CBC (OBC along the 
direction shown with invariance under reflexions assumed), 
J=0.2, and t'=t"=0.0 (600 states, 8 sweeps), (b) Spin-spin 
correlations for 2 mobile holes projected at their most prob- 
able relative distance (circles) in the 2 holes ORBA ground 
state of a 4x6 cluster, J=0.2, t'=-0.35, and t"=0.25. Lines 
indicate AF correlations (thickness proportional to absolute 
value), (c) Same couplings, cluster, technique, and conven- 
tions as in (b) but using PBC (OBC) along legs (rungs), and 
3 holes. Also shown are the hole density-density correlations 
along a center leg, showing that there is no charge order, (d) 
AF spin correlations for a large weight ground-state configu- 
ration of an exactly solved 3 hole, 18 site PBC cluster, J=0.4, 
t'=-0.20 and t"=0.14. 

The half-doped stripes reported here also lead to spin 
IC. For example, in Fig.le the spin structure factor is 
shown for the cases of Figs. la, b. The peak deviation 
from (7r,7r) appears in a robust region of parameter space. 
The spin IC is understood calculating spin-spin correla- 
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tions when two holes in, e.g., the cluster of Fig.lc are 
projected into their most probable location (Fig. 2b): a 
TT-shift across-the-stripe can be clearly observed. The 
across-the-stripe AF correlation strength increases reduc- 
ing J/t and/or increasing t'<0 and t">0 in magnitude. 

Results compatible with n/i^O.5 stripes and associated 
TT-shifts appear in other clusters as well. On a cylin- 
drical 6x4 cluster with PBC along the long direction, 
the 3-holes ground-state has characteristics compatible 
with a doped one-dimensional (ID) closed loop along the 
PBC direction, with 7r-shifts across-the-stripe (see Fig. 2c 
where one of the two degenerate most dominant ground- 
state hole configurations is shown). A h-s-h-s-h-s loop 
(h=hole, s=spin) provides a pictorial representation of 
our results, but this configuration is not rigid neither 
along nor perpendicular to the loop. Density correlations 
along the stripe (Fig. 2c) are actually compatible with a 
ID nft^O.5 system at large on-site U interactions [ p^ , 
suggesting that the stripes described here are metallic. 
No indications of a charge-density-wave along the stripe 
were found. Note also that spin IC induced by antifer- 
romagnetism across holes also exist along the stripes, 
with wavevector 7r/2 for a half-doped stripe. This spin 
IC appears also in half-doped ID models ||l9| . For an iso- 
lated CuO plane, IC should be present in both directions, 
although with quite different wavevectors and intensities. 

Similar results are found in small square clusters: in 
the 2-holes 4x4 lattice with CBC, a 2-hole stripe forms 
along the PBC direction ||]. With PBC in both direc- 
tions, the ground-state resembles a mixture of stripes 
along both axes and since nonzero t'-t" avoids PS, our 
results are not expected to have the boundary effects re- 
cently discussed Tendency to stripe formation is 
found even in tilted clusters: the PBC -s/TSx-s/TS lat- 
tice allows for n/j~0.5 closed loops with 3 holes and such 
structure has a large ground-state weight (Fig. 2d) [ pO| . 
Precursors of the spin structures in Figs. 1,2 appear on 
2- and 3-leg ladders as well, e.g. in Fig. 3a the 2 holes 
ground-state dominant hole configuration of a 3x6 clus- 
ter is shown, with its spin correlations. On 2-leg ladders 
with many holes, 7r-shifts appear at small J/t (Fig. 3b), 
and each hole is "confined" to a rung, precursor of a rung 
stripe as the leg number grows. Spin IC is here found 
both for the 2-leg (Fig. 3c) and 3-leg ladders. 

The results thus far suggest that doped stripes can 
form in spin and hole models using realistic couplings. 
To gain insight into the mechanism driving this complex 
structure, consider now the one hole problem. Fig. 3d 
shows 4-leg ladder spin correlations around a mobile hole 
for momentum (7r,7r). The AF correlations across-the- 
hole are clearly similar to the correlations around the 
individual holes composing the stripes. The 7r-shift char- 
acteristic of the stripes exists in the one-hole state not 
only at (7r,7r) but at several momenta, and, in this sense, 
the spin IC exists already at the one-hole level, a re- 
markable result. Similar conclusions are reached for 3- 
and 2-leg ladders (Fig.3e), and other momenta such as 



(0,7r). Also on small square clusters robust across-the- 
hole AF correlations exist for one hole. Although spin 
IC was found in early studies of the t-J model Q, and 
the nontrivial structures as in Fig. 3d were noticed before 
Q, it was only recently tentatively explained |^ as (lo- 
cal) spin-charge separation, similar to the ID Hubbard 
model where spins across holes are antiparallel [ p^ . 
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FIG. 3. (a) AF spin correlations for 2 holes on the 3x6 
cluster solved exactly at J=0.2, t'=-0.35, t"=0.25, with holes 
projected at their most probable distance in the ground state, 
(b) AF spin correlations at the center of a 2x32 cluster 
with 12 holes (OBC-legs) using DMRG at J=0.2, t'=t"=0.0. 
Shaded regions contain the holes, (c) S(q2:,7r) vs q^; for 
a DMRG 2x16 cluster with 12 and 4 holes and couplings 
(J=0.4, t'=-0.35, t"=0.25) and (J=0.2, t'=t"=0.0), respec- 
tively, (d) Exact AF spin correlations of the PBC 4x6 clus- 
ter with one hole and q=(7r,7r), J=0.2, t'=-0.35, and t"=0.25. 
An elementary block conjectured to form part of stripes is 
framed, (e) Results for 1 hole as in (d) but for a 2x12 clus- 
ter at J=0.2, and a 3x6 cluster at J=0.1, both for t'=-0.35, 
t"=0.25 and q=(0,7r). 

The results shown here lead us to believe that the ob- 
served doped stripes are made out of one-hole building- 
blocks (Fig. 3d). In this respect the insulator limit al- 
ready carries the essential information needed to build 
the stripes, providing an unexpected potential simple link 
between undoped and doped cuprates. This is compati- 
ble with the behavior of the large energy scale PES pseu- 
dogap which can be traced back to the one-hole disper- 
sion of the insulator |ll|, suggesting a smooth evolution 
from the undoped to underdoped regimes. 

However, further elaboration is needed since for one- 
hole the lowest energy is found at q~(7r/2, 7r/2) ]lO| , p^ . 
Naively, hole pockets at (7r/2,7r/2) should appear at fi- 
nite hole density. In addition, across-the-hole AF bonds 
are weaker at (7r/2,7r/2) than at momenta such as (7r,0) 
or (7r,7r) although they are still present. To address 
this issue let us calculate (nq) = (c]jCq), i.e. the ground- 
state hole number with a given momentum q (note that 
(uq) includes both coherent and incoherent weight). As 
example, consider the two- hole problem on the 4x6 lat- 
tice of Fig. 2b. The interesting result in Fig. 4a is that 
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FIG. 4. (a) (nq) vs q, for the 2 hole ground-state of the 
4x6 cluster of Fig. 2b. (b) Qualitative representation of a 
one-hole state with strong AF correlations across-the-hole as 
the building-block of nh~0.5 stripes. In the 1 hole case the 
frustration effect is shown with question marks, (c) Exact 
spectral function A(q, oj) for one-hole on the 4x6 cluster with 
PBC in both directions, J=0.2, q=(7r,7r) and t7t"=-1.4. Val- 
ues of t' as well as the (small) weight in the first pole (and 
its location) are indicated. Note the accumulation of weight 
at large energies, (d) Same as (c) but for a 2x12 cluster with 
PBC along legs and q=(0,7r); (e) Same as (c) but for a 3x6 
cluster (PBC-leg, OBC-rung) at q=(0,7r). 

the ground state carries dominant weight at momenta 
around (tTjTt), and the one- hole states with this mo- 
mentum have robust AF correlations across-the-hole 
(Fig. 3d), compatible with our conjecture There 
are no indications of small hole-pockets in our studies, 
and the Fermi surface appears open. In this frame- 
work the across-the-hole correlations of the, e.g., (7r,7r) 
or (7r,0) holes can be "linked", as pictorially shown in 
Fig. 4b, improving the hole mobility since now they share 
a large region where they do not need to fight against 
the spin background to move. Creating a stripe loop 
also avoids the spin frustration intrinsic of the individ- 
ual hole states when across-the-hole robust correlations 
are present (Fig. 4b). In addition, our results help under- 
standing better the observed stripe density: for n/i^l the 
across-the-hole AF bonds in the stripe direction cannot 
form and holes do not improve their kinetic energy, while 
for a very hole diluted stripe the finite-size elementary 
blocks (Fig. 3d) do not touch and cannot have a com- 
mon spin arrangement. For completeness, in Figs. 4c, d 
the one-hole spectral function is exactly calculated on 
4-, 3- and 2-leg ladders at small J/t. Note the remark- 
able small quasiparticle weight, correlated with a robust 
across-the-hole AF correlation (see also |^). The one- 
hole states contributing to stripes have exotic properties, 
including a tendency to spin-charge separation |^ . 

Summarizing, indications of n/j~0.5 stripes were found 
in the extended t-J, t-J^, and (at small J/t) in the stan- 
dard t-J models. The gain of kinetic energy against the 



loss of AF energy appears enough to stabilize stripes, 
namely the driving force is a one hole process and the 
seed for stripes is already present in the insulator. Con- 
trary to most approaches to stripe formation, here the 
small J/t regime was emphasized. The scenario reported 
here is a generalization of the ID spin-charge separation 
involving individual holons, with the twist that stripes of 
holons are needed in 2D to avoid frustration. This result 
is compatible with Zaanen's picture of stripes as "holons 
in a row" Charge and spin could be separated in 2D 
in more subtle ways than anticipated. 
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